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Abstract 

Jet production in deep inelastic scattering for 120 < Q 2 < 3600 GeV 2 has been stud- 
ied using data from an integrated luminosity of 3.2 pb _1 collected with the ZEUS de- 
tector at HERA. Jets are identified with the JADE algorithm. A cut on the angular 
distribution of parton emission in the 7*-parton centre-of-mass system minimises the 
experimental and theoretical uncertainties in the determination of the jet rates. The 
jet rates, when compared to 0(a s 2 ) perturbative QCD calculations, allow a precise 
determination of a s (Q) in three (^-intervals. The values are consistent with a run- 
ning of a s (Q), as expected from QCD. Extrapolating to Q = M z o yields a s (M z o) = 
0.117 ± 0.005 (stat) ±8:Soi (syst exp ) ± 0.007 (syst theory ). 
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1 Introduction 



Neutral current (NC) deep inelastic scattering (DIS) (Ip — ► IX; I — e, fi), is characterised by the 
exchange of a virtual photon or Z° boson between the incident lepton and proton. In the naive 
quark-parton-model (QPM) the process V*q — > q (V = 7, Z°) gives rise to 1+1 jets in the final 
state corresponding to the struck quark from the proton and the proton remnant (hereafter 
denoted by "+1"). Multi-jet production in DIS beyond 1+1 jets provides a good laboratory 
for testing quantum chromodynamics (QCD). From the measured rate of 2+1 jet events it is 
possible to determine the strong coupling constant a s , for fixed kinematics and a given jet 
definition, by comparing to theoretical calculations in which a s is the only free parameter. 

To leading order in a s , 2+1 jet production proceeds via QCD-Compton scattering (V*q — > qg) 
and boson-gluon fusion (BGF) (V*g — > qq). For the extraction of a s from the measured 
jet rates to be reliable the 2+1 jet rate must be calculated at least to next-to-leading order 
(NLO) in QCD, where the renormalisation scheme is defined unambiguously. Furthermore the 
jet definition has to be treated in the same way in theory and experiment for a quantitative 
comparison with the predictions of QCD. Theoretical calculations |], @, f|| for the jet rates at 
the parton level are currently available only for the JADE jet definition scheme |>J. Therefore 
the measured jet rates, obtained using the same jet-finding scheme, have to be corrected to 
the parton level so that a comparison with the NLO 0(a s 2 ) calculations can be made in 
order to determine a s . The extracted a s value can be expected to be reliable when the NLO 
calculations reproduce the corrected jet rates over a wide kinematic range and the extracted 
value is insensitive to the cuts applied at the detector level. In this analysis a cut on the parton 
variable z (described later) is applied, which restricts the phase space so that these requirements 
are well satisfied. 

Multi-jet production in DIS has been studied by the E665 fixed-target experiment at FERMI- 
LAB at a low centre-of-mass energy, ^/S, of ~ 30 GeV ||, and at higher energies, -y/i=300 GeV, 
by ZEUS 0] and HI || at HERA where jet structures are more clearly discernible. This pa- 
per describes the extraction of a s from measurements of multi-jet rates at Q 2 between 120 
and 3600 GeV 2 . An earlier study of jet rates and jet kinematics has been reported by this 
experiment ||; an extraction of a s from multi-jet production has been reported by HI ||10| . 



2 The ZEUS Detector 

The data used in this analysis were collected with the ZEUS detector during 1994 when HERA 
provided collisions between 27.5 GeV electrons or positrons^ and 820 GeV protons, yielding a 
centre-of-mass energy of 300 GeV. They correspond to an integrated luminosity of 3.2 pb _1 . 

ZEUS is an almost hermetic, multipurpose, magnetic detector and has been described elsewhere 
in detail ||11|| . Here a brief description of the components relevant for this analysis is given. 
Charged particles are tracked by the inner tracking detectors which operate in a magnetic field 
of 1.43 T provided by a thin superconducting coil. Immediately surrounding the beam pipe 
is the vertex detector, a drift chamber, which consists of 120 radial cells, each with 12 sense 
wires ||12|| . It is surrounded by the central tracking detector which consists of 72 cylindrical 



1 Hereafter "electron" is used in a generic sense to refer to e 
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drift chamber layers, organised into 9 'superlayers' [13]. In the present analysis these tracking 
detectors are primarily used for the determination of the event vertex. 



The energy associated with the hadronic final state and the scattered electron is measured 
with the uranium-scintillator calorimeter (CAL) [|14]] which consists of three parts: the forward 
(FCAL), the rear (RCAL) and the barrel calorimeter (BCAL).The ZEUS coordinate system is 
defined as right handed with the Z axis pointing in the proton beam direction, hereafter referred 
to as "forward" . The X axis points horizontally towards the centre of HERA and the Y axis 
points vertically upwards. The polar angle 9 is defined with respect to the Z direction. Each 
part of the calorimeter is subdivided longitudinally into one electromagnetic section (EMC) and 
one hadronic section (HAC) for the RCAL and two HAC sections for BCAL and FCAL. Holes 
of 20 x 20 cm 2 at the centre of FCAL and RCAL accommodate the HERA beam pipe. In the 
XY plane around the FCAL beam pipe, the HAC section is segmented in 20 x 20 cm 2 cells and 
the EMC section in 5 x 20 cm 2 cells. In total, the calorimeter consists of approximately 6000 
cells. In terms of pseudorapidity, r] = — In tan |, the FCAL covers the interval 4.3 > i] > 1.1, 
the BCAL 1.1 > r\ > —0.75 and the RCAL —0.75 > r] > —3.8, for the nominal interaction 
point at X = Y = Z = 0. The CAL energy resolution, as measured under test beam conditions, 
is oe/E = Q.l&/\/E~ for electrons and (Je/E = 0.35/v^ for hadrons (E in GeV). The time 
resolution of the calorimeter, which is important for rejecting beam-gas backgrounds, is better 
than 1 ns for energy deposits greater than 4.5 GeV. 



3 Event Kinematics 

For a given ep centre-of-mass energy y^s, the differential cross section for leading order 0{a s 1 ) 
2+1 jet production in DIS depends on 5 independent kinematic variables, which we take as 
x, y, x p , z, and $ Hl5fl . The first two, Bjorken-x and y, are sufficient to describe the O(a s ) 
QPM 1+1 jet process. They correspond to the momentum fraction of the proton carried by 
the struck quark (x) and the fractional energy transfer between the electron and the proton in 
the proton rest frame (y). Three additional variables (x p , z, $) are introduced to describe the 
2+1 parton kinematics. The parton variable x p is defined by: 

Q 2 Q 2 x 

Xp 2p-q Q 2 + m% £' 

where q is the four-momentum of the exchanged virtual boson in the ep scattering process, £ 
is the fraction of the proton's four-momentum P carried by the incoming parton with four- 
momentum p = £P, rriij is the invariant mass of the two non-remnant jets and Q 2 = — q 2 . Q 2 , x 
and y are related by Q 2 = s x y. The parton variable z is defined by: 

z\ = = -■(! — cos 6^, 



p-q 2 v l/ Ei=i, 2 cos 0*) 

The formula is given for one of the partons i = 1. The outgoing four-momentum of the parton 
from the hard scattering is p\ and d\ is the scattering angle in the 7*-parton centre-of-mass 
system. Experimentally, z is determined in the HERA system from the energies and angles, 
Ei and 0j, of the two jets. The jets are assumed to be massless. The other parton satisfies the 
constraint z 2 = 1 — Z\. The angle $ represents the azimuthal angle between the parton and 
lepton scattering planes in the 7*-parton centre-of-mass system. 
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Since the ZEUS detector is nearly hermetic, it is possible to reconstruct the kinematic variables 
x, y and Q 2 for NC DIS using different combinations of the angles and energies of the scattered 
lepton and of the hadronic system (T^|. The electron method was used to determine y as 
y e from E' e and 6 e , the energy and polar angle of the scattered electron. The hadronic, or 



Jacquet-Blondel method flT7| , was used to reconstruct y as yjs = E/i-^/i(l — cos8h)/(2E e ) 
where and 9^ are the energy and polar angle calculated from the calorimeter cells not 
associated with the scattered electron, and E e is the electron beam energy The double angle 
(DA) method uses 6 e and 7^, the polar angle of the struck quark in the QPM which is given 
by cos7# = {T,hPr,h ~ ( 2E eyjB) 2 ) /(EhPT,h + (2E e y JB ) 2 ). The DA method, which measures 
Q 2 with small bias and good resolution in the kinematic range of this analysis, was used to 
reconstruct the x, y and Q 2 variables of NC events (and the jet variables x p and z defined 
above) fllT 



4 Trigger Conditions and Event Selection 

The trigger and event selection followed closely that described in reference |J. The trigger 
acceptance was essentially independent of the hadronic final state with an acceptance greater 
than 97 % for Q 2 > 10 GeV 2 . Neutral current DIS events were selected by the following 
criteria: the event times measured by the FCAL and the RCAL had to be consistent with 
an interaction inside the detector. This cut strongly reduced beam-gas background. The Z 
position of the event vertex was reconstructed from the tracking data. Events were accepted if 
the Z position was within ± 50 cm of the nominal interaction point. An electron candidate with 
energy greater than 10 GeV had to be found in the calorimeter. To reject backgrounds from 
photoproduction events with a fake electron (mostly 7r°'s close to the proton beam) the electron 
candidate was required to satisfy y e < 0.95. Photoproduction and beam-gas backgrounds were 
further suppressed by demanding energy- momentum conservation. For a fully contained event, 
and neglecting the detector resolution, one expects E — P z = 2 ■ E e , where E and Pz are the 
summed energy and Z-component of the momenta of all objects measured in the calorimeter. 
Taking the detector resolution into account 35 GeV< E — P z < 60 GeV was required to select 
DIS events. The background from photoproduction was estimated to be negligible and that 
from QED Compton scattering was found to be less than 1%. Diffractive events, which do not 
deposit a significant amount of energy in the FCAL, did not pass the selection criteria given 
below. Finally, beam halo muons and cosmic rays were rejected by suitable algorithms. 

Several considerations motivated the selection of the kinematic region used for the a s measure- 
ment. First, the analysis was restricted to high Q 2 , where clear jet structures are observed 
and hadronisation uncertainties are minimised. Secondly, theoretical uncertainties in the 2+1 
jet cross section are small at high x, where the parton densities of the proton are well known. 
In addition at high x the uncertainty stemming from the initial state parton-showers, used in 
the Monte Carlo simulation to correct the data to the partonic level, was reduced. Thirdly, 
the acceptance for 2+1 jet events increases at high y: in particular, the forward jet is well 
contained within the detector. Finally, all of the above concerns were balanced against the 
need for sufficient statistics. The kinematic region selected for the final analysis was therefore: 
120 < Q 2 < 3600 GeV 2 , 0.01 < x < 0.1, and 0.1 < y < 0.95, resulting in a sample of 4472 
events. The Q 2 range was further subdivided into three regions to measure a s (Q) at increasing 
scales as a consistency check and as a test for the running of the strong coupling constant. 



3 



These ranges were: 120 < Q 2 < 240, 240 < Q 2 < 720, and 720 < Q 2 < 3600 GeV 2 . The 
number of events in each region were 1649, 2048 and 775, respectively. 



5 Jet Definition and Jet Kinematics 

The JADE algorithm |J was used to relate the hadronic final state measured in the detector to 
the underlying hard scattering processes. It is a cluster algorithm based on the scaled invariant 
mass, i/ij = m 2 j/W 2 = 2EiEj(l — cos8ij)/W 2 , where m^- is the invariant mass of the two 
objects i and j, which are assumed to be massless. The scale W 2 is the squared invariant 
mass of the overall hadronic final state and E{,Ej and 9ij are the energies of the objects and 
the angle between them. Starting with the minimum of all possible combinations, objects 
were merged by adding their four-momenta until for all objects exceeded a jet resolution 
parameter y cut . Those objects remaining were then considered as jets. The JADE algorithm 
was slightly modified |18|, [19| for use at the detector level in ep collisions by the addition of a 
pseudo-particle inserted along the Z axis. The missing longitudinal momentum in each event 
was assigned to the momentum of the pseudo-particle. It prevents the detected fraction of 
particles originating from the proton remnant from forming spurious jets. 

The measured calorimeter energies above 150 (200) MeV for EMC (HAC) cells and their angles 
relative to the interaction point were used to define vectors which were input to the JADE 
algorithm in the detector level analysis. At the detector level, the scale W was calculated as 
Wvis = s (1 ~ x da) Ujb- This value reflects the measured rather than the true hadronic activity 
and so reduces the event-by-event correction for the detector resolution when calculating yij. 

For 2+1 jet production in DIS, one of the non-remnant jets is typically directed forward because 
of the forward singularity in the cross section. Such forward singularities are regulated in the 
theory by the cutoff y cut . Requiring a large y cut is, however, not sufficient to avoid the problems 
arising from forward-going jets close to the beam pipe and the proton remnant. This is achieved 
by a cut on the parton variable, z, for 2+1 jet events^. In QCD the rapid rise towards z = 
results from collinear and infrared singularities. In order to avoid this kinematical region the 
analysis was restricted to events satisfying 0.1 < z < 0.9. This requirement also reduces the 
fraction of forward jets (9j et < 8°) from 30% to 10%. Figure |l|a shows the dR 2 +i/dz distribution 
for 2+1 jet events. Here Rj+i = Nj + i/N to t, where j stands for 0, 1, 2, or 3, Nj + x is the number 
of (j+1) jet events and N to t is the total number of selected DIS events. Figures |l]b d show the 
resulting x p , px and distributions. The predictions of the NLO calculations of DISJET and 
PROJET (discussed later) are also shown in Fig. |l| The z-cut results in a significantly improved 
agreement between the calculations and the data compared to our earlier analysis done without 
this restriction on z 0. This cut removes jets with transverse momenta pt below ~ 4 GeV 
where pr is measured with respect to the 7* direction and is calculated in the 7*-parton system 
as: 



2 z is not denned for 1+1 and 3+1 jet events. In the following, 2+1 jet events failing to pass the z cut are not 
considered as 2+1 jet events. 1+1 and 3+1 jet events are counted in Ntot-, the total number of selected events. 
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6 Jet Reconstruction and Jet Rates 



The acceptance and resolution, as well as the correction of the measured jet rates to the parton 
level, were determined using Monte Carlo methods. Neutral current DIS events, generated 
using LEPTO 6.1 [^] and the Lund string fragmentation model |2l| for the hadronisation, 
were interfaced via DJANG06 2.1 p[ to HERACLES 4.1 J231 for QED radiative corrections. 



They were passed through a GEANT |24| based detector simulation, and subsequently analysed 
with the same reconstruction, selection and jet analysis procedures as the data. Both the hard 
emission of partons at the matrix element level (calculated to leading order in a s ) and the higher 
order soft parton showers are included in the LEPTO matrix element, parton shower (MEPS) 
model. The MEPS model satisfactorily describes the global jet properties and production rates 
observed for the data in the selected kinematic region Q. 

When generating events with the MEPS model, default values of all parameters were used 
except for the parameter y m i n , which sets a minimum y^ of partons in first order QCD matrix 
elements [20|. The value of y m in was lowered from 0.015 to 0.005 in order to study the measured 
jet rate as a function of the jet resolution parameter y cut for y cut >0.01. The parton densities 



of the proton were taken from the MRSD'- set ||25|| . 

The jets were reconstructed by applying the JADE algorithm at the parton level, the hadron 
level and the detector level. These jets were constructed respectively from the output of the 
parton shower step of the event generator, the true momenta of the hadrons before the detector 
simulation and the energy deposits in the calorimeter cells after the detector simulation. The 
ratio of the jet rates at the different levels of the event simulation are the corresponding cor- 
rection factors for hadronisation (Ch) and detector simulation (Co), with which the measured 
jet rates were multiplied. Both the detector and hadronisation corrections were found to be 
below 20%. Table [3] shows the correction factors and the corrected jet rates -R2+1 for the three 
Q 2 intervals and for the combined Q 2 region. 
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1.07 


3.5±0.3 


0.99 


1.10 


0.060 


2.3±0.4 
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Table 1: 2+1 jet production rates (in %) corrected to the parton level (R2+1) and correction 
factors for detector effects (Cd) and for hadronisation (Ch) in the three Q 2 intervals and for the 
combined region. Errors shown are statistical only. 



7 0(a s 2 ) Perturbative QCD Calculations 

For NC electron-proton scattering, the 2+1 jet differential cross section at the 0(a s 1 ) LO level, 
expressed in terms of the above variables, is given by: 
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dx dy y Q 2 J x p J J 2n 



where I g and I q are the gluon- and quark-initiated contributions respectively |15| which contain 
singularities at z — 0, z = 1 and x p — 1. The singularities correspond to the limit where two 
partons are irresolvable due to vanishing energy of one of the partons or vanishing opening 
angle. In the BGF process the singularity is related to the collinear emission of an outgoing 
quark in the proton remnant direction; for the QCD Compton process the singularity occurs 
when the momentum of the gluon is parallel to that of the quark or when the quark emits a 
very soft gluon. The integration must be done numerically because of the x-dependent parton 
density distributions. The integration limits for z and x p in the JADE scheme are functions of 
the scaled invariant mass cutoff y cu t = m 2 j/W 2 , where W is the reference mass scale and rriij 
is the invariant mass between any two partons [F2J. Any pair of partons with a scaled invariant 
mass below this cutoff is not resolved. Therefore the singularities are regulated by a single 
cutoff parameter y cut . The leading order (LO) and next-to-leading (NLO) order 2+1 jet cross 
section can be expressed as: 

a<J 2+l 

- — — = c 3 i • a s , 

dx dy 

d(T NLO 

2+1 2 

■_- = Csi .a. + c B .a.. 

Next-to- leading order corrections to d<72+i include the contribution from unresolved 3+1 jet 
events as well as negative corrections coming from virtual loops ||, ||. The coefficients c^- 
contain the hard scattering matrix elements and the parton density functions of the incoming 
proton. The effect of a change in a s on the parton densities is negligible for our Q 2 range [26 . 
The first index, i, stands for the jet multiplicity (including the remnant jet) and the second 
index, j, represents the order of the a s calculation. After integrating over the jet variables 
(x p , z, $) the coefficients c^- are functions of the event kinematic variables x, y, y cut and the 
factorisation scale fip; c 32 depends also on the renormalisation scale fiR pTJ. The parton 



densities contained in are calculated at the scale Q 2 . In finite order perturbative QCD 
calculations a s depends on the renormalisation scale /xr. The 2+1 jet rates are derived from 
the cross sections by R2+1 = (?2+i/<J to t- The resulting 0(a s 2 ) corrections to -R2+1, using the 
NLO calculations, are considerable: they vary from -20 to +20% when y cut is varied between 
0.01 and 0.06 in the kinematic region used for this study [§], ||. The numerical cross section 
calculations are available in the DISJET program |28| by Brodkorb and Mirkes, and in the 
PRO JET program |29| by Graudenz. Both programs agree in their predictions of a s for a given 
jet rate and they reproduce the shape of the measured jet rate distributions as a function of 
Ucut well in the investigated kinematic range (see below). 



The renormalisation scheme used in the calculation is the MS scheme. In second order, the 
dependence on other renormalisation schemes can be completely specified by one parameter, 
which can be chosen to be the value of the renormalisation scale /Jr. We chose ijl 2 r = Q 2 for 
our analysis. The same scale is chosen for the factorisation scale fiF- The parton densities 
were calculated with a fixed Ap^- = 154 MeV. In the kinematic range used in this analysis, the 



effect of varying A^- in the parton densities is expected to be small [p6 |. In the programs the 
contributions from the c— and b— quarks are zero in the parton density parametrisations below 
the single quark mass thresholds, as defined in the MS factorisation scheme. Above threshold, 
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the contributions from the c— and b— quarks are calculated assuming zero mass. The number 
of flavours in the formula for the running coupling constant is changed at the single mass 
threshold, as required by the MS renormalisation scheme, giving rise to five flavours for Ajt^ 
if Q 2 > m 2 . At the BGF vertex we have used five flavours too because in the kinematic range 
of this analysis mfj is above 4 • m 2 . Using four flavours at the BGF vertex in the PRO JET 
program would increase the a s (M z o) value by 0.0025. In the x, Q 2 region under study the 
contribution from massive b— quarks to the proton structure function is calculated to be below 
2% pll. 



8 Determination and Q 2 Dependence of a s 

The value of a s was determined by varying in the QCD calculation until the best fit to the 
ratio R2+1 was obtained at y cut = 0.02. The slope of the measured -R2+1 as a function of y cut 
agrees with that from the calculation, showing that the result is not sensitive to the particular 
value of y cut used. We chose y cut = 0.02 for the fit because the contribution from R3+1, which 
is a higher order effect, becomes negligible for y cut > 0.02. Furthermore the statistics are large 
and the jets are resolvable at this value and the 2-jet system has a large invariant mass. 

Figures 0a-d show the corrected jet rates, i?i+i, R2+1 (also shown in Tab. [I]) and -R3+1 as a 
function of y cut for data compared with the DISJET and PROJET NLO QCD calculations for 
the three Q 2 intervals and for the combined region. Only statistical errors are shown. All NLO 
terms are taken into account in both programs; however, they use different approximations 
for some of these terms. There is good agreement between the corrected jet rate and the 
NLO QCD calculation over most of the range in y cut shown and in particular at the nominal 
Vcut — 0.02, where the a s value was extracted for this analysis. Both programs agree well 
in their prediction of the jet-rate dependence as a function of y cut . The best fit values for 
a s are used in the calculation. The range in y cut was restricted to 0.01 to 0.06 because at 
lower values of y cut the jets are not experimentally resolvable and higher order corrections 
are significant, while at larger values terms proportional to y cut , neglected in the calculation, 
become significant. Moreover, uncertainties in the renormalisation scale and the hadronisation 
corrections also become large for y cut above ss 0.06 |27jj . 

The values of a s are plotted in Fig. |3] as a function of Q for each of the three Q 2 ranges. They 
are calculated from the fitted values of A^-. In Tab. the a s values determined for the three 

Mo u 

ranges in Q as well as for the full kinematic range are listed. Also shown are the values of a s 
extrapolated to Q — M z o. Both statistical and systematic errors (discussed in the following 
section) are given. 

(5) 

In addition Fig. shows the curves for A-^ = 100, 200, and 300 MeV. The measured a s 
decreases with increasing Q, consistent with the running of the strong coupling constant if Q 
is taken as the scale. The fit to a running a s (where a s was determined in the full Q 2 range) 
yields a x 2 of 2.2 for 2 degrees of freedom, which corresponds to a confidence level of 58.6%. 
A least squares fit to the hypothesis of a constant a s was performed. Only statistical errors 
were considered in this fit as the systematic errors are strongly correlated. This fit yields a x 2 
of 7.7 for 2 degrees of freedom, which corresponds to a confidence level of 2.1%. Taking into 
consideration the systematic uncertainties, the x 2 f° r constant a s varies from 4.4 (changing 
energy scale by -5%) to 10.3. A constant a s is thus ruled out at 90% confidence level. The 
three values of a s , expressed at the mass of the Z° boson, are consistent within the errors. 



7 



Q 2 

( GeV 2 ) 


< Q > 
(GeV) 


MS 

( MeV ) 


a s (Q) 


a s (M z o) 


19fl 9Afl 


1 3 3 


251 


+108 +31 
-97 -74 


+115 
-105 


0.171 


+0.015 +0.005 
-0.017 -0.012 


+0.016 
-0.018 


0.120 


+0.007 +0.002 
-0.008 -0.006 


+0.007 
-0.009 


240 - 720 


20.4 


217 


+90 +76 
-74 -60 
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-0.006 
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35.5 


86 
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-58 -47 


+61 
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+0.013 +0.006 
-0.017 -0.012 


+0.010 
-0.006 


0.103 


+0.010 +0.004 
-0.013 -0.010 


+0.008 
-0.004 


120 - 3600 


22.1 


208 


+64 +57 
-53 -50 


+89 
-75 


0.148 


+0.008 +0.007 
-0.008 -0.007 


+0.011 
-0.012 


0.117 


+0.005 +0.004 
-0.005 -0.005 


+0.007 
-0.007 



Table 2: The measured values of and a s for the three ranges in Q 2 as well as for the full 
Q 2 range. The first error is statistical, the second corresponds to the experimental systematic 
uncertainty and the third to the theoretical systematic uncertainties (hadronisation, parton 
density and scale uncertainty). 



9 Systematic Uncertainties 



9.1 Experimental, Hadronisation and Parton Density Effects 

Sources of systematic uncertainties in the a s determination were grouped into the following 
classes: event selection, energy scale, jet analysis, fitting method, model dependence of detector 
corrections, hadronisation corrections and parton density (see Fig. |4]). The first five classes were 
attributed to the experimental systematic error. The uncertainties were studied for each of the 
three Q 2 ranges separately as well as for the combined kinematic region. Only the systematic 
uncertainty from the latter study is described here. To illustrate the systematic uncertainty 
in a s associated with each systematic effect, the fitted a s value obtained when the systematic 
effect is varied was compared with the central value of a s (see Tab. |2|). 

The subgroups of experimental systematic errors are denoted by (a)-(e). The systematic errors 
from the event selection (a) included: effect of using a different electron finding algorithm; 
variations of the selection criteria, E — Pz > 45 GeV; y e < 0.7. The errors from the energy 
scale (b) included a ±5% error assigned to the uncertainty of the calorimeter energy response. 
The errors from the jet analysis (c) included: the choice of a different mass scale, W^ A = 
s (1 — xda) Vda and Wj B — s (1 — xjb) Vjb, in the JADE scaled mass definition, = m\-jW 2 \ 
the cells around the FCAL beam pipe, which contain mainly the proton remnant, were first 
preclustered and the resulting objects were used in the jet clustering algorithm (instead of the 
cell vectors themselves). The errors from the fitting method (d) included: a QCD fit at y cu t=0.03 
instead of y cut =0.02; the analysis was cross-checked by a QCD fit to the differential jet rates, 
D 1+1 , defined by D 1+1 (y cut ) = [R 1+1 (y cut + Ay cut) - R 1+1 (y cut )]/ Ay cut ; a more restrictive z-cut, 
0.15 < z < 0.85, was used. Finally, the error from the model dependence of the corrections for 



the detector acceptance and resolution (e) was estimated by using the colour-dipole model [[H 
as implemented in the ARIADNE 4.06 Monte Carlo ||32|| . The largest uncertainties for each 
subgroup were added in quadrature to give the experimental systematic error. 

To evaluate the uncertainty of the hadronisation correction, several aspects of the hadronisation 
scheme were varied while the standard detector corrections based on the LEPTO MEPS Monte 
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Carlo were retained. These studies were performed at the generator level. First, parameters in 
the Lund string fragmentation model |$3| were varied: a in the 'symmetric fragmentation func- 



tion', which regulates the longitudinal quark fragmentation, was varied between 0.1 and 1; <jp t , 
which controls the hadron transverse momentum distribution was varied between 0.25 and 0.45 
GeV. Second, parameters of the parton shower model employed in the LEPTO MEPS Monte 
Carlo were changed: y m in was varied between 0.005 and 0.015; the minimum virtuality scale, 
Qo, at which the parton showering is stopped, was changed from 0.8 to 4 GeV; the primordial 
transverse momentum h? of the struck parton in the proton was varied from 0.44 to 0.7 GeV. 
Finally, a completely different hadronisation model as implemented in the HERWIG 5.8 Monte 
Carlo p4| was used. Most of these changes result in relatively small systematic errors in a s as 



shown in Fig. f|. The two largest deviations from the central value of a s arise from the change 
of the hadronisation model and from the variation of Qq. 

We also repeated the analysis with parton density sets MRSA, GRV HO, and CTEQ 3M in 
the NLO calculation, all of which describe the results from present DIS data well [16[]. The 



differences in a s (22.1 GeV) from the central value are small (< 0.0022), as shown in Fig. The 
fitted a s value depends only weakly on the a s value used in the parton density parametrisations 



In the x, Q 2 region under study the contribution from massive b— quarks to the proton structure 



function is calculated to be below 2% [3(J. The effect of calculating with four instead of five 
flavours at the BGF vertex was estimated with the PRO JET program and was found to increase 
a s (Mzo) by 0.0025. This number is not included in the systematic errors given. 



9.2 Scale Dependence Effects 

Our best estimate of the scale uncertainty in the measured a s was obtained from DISJET and 
PROJET by varying //r and \ip from 0.4 Q 2 to 2.0 Q 2 , redoing the fit to the jet rates and 
evolving to obtain the corresponding value of a s at the original scale Q (shown in Fig. ^ for 
the full Q 2 range). The scale dependence decreases with increasing Q and becomes negligible 
in the highest Q 2 interval. It is slightly larger in DISJET than in PROJET. 

Deep inelastic production of jets is a multi-scale process and it is not evident that Q 2 is the best 
choice 0, |35| for the renormalisation and factorisation scales in the perturbative calculation. 
Alternative scales have been suggested, e.g. p\ of the jets or the square of the invariant mass, 
m 2 j, of the two jets. As a simple test the ratios < p%/Q 2 > and < m 2 -/Q 2 > were evaluated 
for the full Q 2 range for 2+1 jet events and were used to estimate the resultant change of scale 
and hence of the uncertainty in a s . For our events these ratios typically lie between 0.4 and 2, 
i.e. within the range explored above in our estimation of the scale uncertainty using DISJET 
and PROJET. 

For each group in Fig. ^ we quote the largest deviation from the central value in each direction 
as the systematic error. The positive and negative deviations are then added in quadrature 
separately to give the systematic error. The total systematic uncertainty for the value of a s 
resulting from the effects studied in Fig. [| is comparable to the statistical errors. In the final 
result given below the uncertainties from the experimental and theoretical systematic effects 
(hadronisation, parton density distributions and scale effects) are quoted separately. 
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10 Summary 



Multi-jet production in ep collisions was investigated using the JADE jet definition. In ep 
collisions the application of the single jet resolution parameter, y cut , is not sufficient to restrict 
the phase space of 2+1 jet production to an experimentally well understood and theoretically 
safe region. An additional cut on the parton variable z is introduced, which excludes the 
problematic region where higher order effects are important and jets are not well measured in 
the experiment. With this additional cut the multi-jet production rate in DIS is well reproduced 
by 0(a 2 ) perturbative QCD calculations. 

The value of a s (Q) was determined in three Q 2 regions in a single experiment and was found 
to decrease with Q, consistent with the running of the strong coupling constant. 

The value for a s using the data from the entire kinematic range 120 < Q 2 < 3600 GeV 2 , and 
expressed at the Z° mass is given by: 

a s (M z o) = 0.117 ± 0.005 (stat) t°;°05 (exp) ljg« (had) ±° ™{ (pd) ±g;ggg (scale) 
= 0.117 ± 0.005 {stat) t° Zt (syst exp ) ±0.007 (syst theory ), 

where stat corresponds to the statistical error and the systematic error components (syst) 
consist of the experimental (exp), hadronisation (had), parton density (pd) and the scale (scale) 
related uncertainties. The overall systematic error is separated into its experimental (exp) and 
theoretical (theory) contribution. 

Our value of a s is consistent with the most recent compilation by the Particle Data Group |36| of 
previous measurements of a s (Mzo) using different methods: 0.112 ±0.005 (DIS), 0.121 ±0.006 
(e + e~ event shape analysis) and 0.124 ± 0.007 (Z° width). The good agreement between 
our value of a s and the results obtained using other methods in different kinematic regimes 
represents a significant test of QCD. 
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Figure 1: (a) Distribution of the parton variable, z, of one of the two non-remnant jets 
in 2+1 jet events in the range 120 < Q 2 < 3600 GeV 2 , compared to the NLO calculations 
(PRO JET and DISJET). The dots with error bars are the measured data. The curves represent 
the theoretical predictions after the application of the cuts. The histograms show the same 
theoretical prediction with the binning of the data, (b) Distribution of x p for 2+1 jet events, 
(c) Transverse momentum distribution for the two jets, (d) Invariant mass distribution 
of the two non-remnant jets. Only events satisfying 0.1 < z < 0.9 were plotted in Figs. b-d. 
All jet rates are evaluated for y cut = 0.02. The data points are corrected to the parton level 
and plotted with their statistical errors only. 
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Figure 2: Jet production rates Rj as a function of the jet resolution parameter y cut for Q 2 in 
the range (a) 120 < Q 2 < 240 GeV 2 , (b) 240 < Q 2 < 720 GeV 2 , (c) 720 < Q 2 < 3600 GeV 2 , and 
(d) 120 < Q 2 < 3600 GeV 2 . Only statistical errors are shown. Two NLO QCD calculations, 
DISJET and PROJET, each with the value of obtained from the fit at y cut =0.02, are also 
shown. 
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Figure 3: Measured values of a s (Q) for three different Q 2 regions. The statistical error 
corresponds to the inner bar and the thin bar reflects the statistical and systematic error added 
in quadrature. Note that the systematic errors are strongly correlated. The dashed curves 
represent a s with = 100, 200, and 300 MeV. 
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Figure 4: Systematic uncertainties in the measured value of a s (and A^-) for Q 2 in the 
range 120 < Q 2 < 3600 GeV 2 expressed as the deviation from the central value for the listed 
alterations in the analysis. Sources of systematic uncertainties are grouped into the four areas: 
experiment, hadronisation correction, parton density, and scale. The experimental uncertainty 
is subdivided into: (a) event selection, (b) energy scale, (c) jet analysis, (d) fitting method, and 
(e) model dependence of the detector correction. 
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